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Abstract 

A coniplete analysis of the consequences of introducing a set of holonomic gauge 
fixing constraints (to fix the dynamics) into a singular Lagrangian is performed. 
It is shown in general that the dynamical system originated from the reduced 
Lagrangian erases all the information regarding the first class constraints of the 
original theory, but retains its second class. It is proved that even though the 
reduced Lagrangian can be singular, it never possesses any gauge freedom. As an 
application, the example of n • A = gauges in electromagnetism is treated in full 
detail. 
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1 Introduction 

Take the Lagrangian for electromagnetism, L = —jF^uF'^^'^. L is singular because its 
Hessian matrix with respect to the velocities is singular. Now consider a generic axial 
gauge n • ^ = 0, with ng 7^ 0, and plug this constraint into the Lagrangian to eliminate 
. We end up with a reduced Lagrangian Lr wich turns out to be regular and that 
can be used to define the propagator. The question is: is this procedure (plugging 
the constraint into the Lagrangian) correct?. Or put it in another way: what is the 
relation between the dynamics defined by L and that by L/j?. Also: what happened to 
the Gauss Law V • E = 0, which is a constraint for L but is no longer present for L_r?. 

In this paper we will give an answer to these questions by studying the relations 
betwwen L and L/j in the general case of a Constrained Dynamical System (CDS). 
There is sometimes a bit of confusion when using these words. Here we mean dynami- 
cal systems defined through a LagrangianQ that happens to be singular (i.e.: when the 
Hessian matrix of the Lagrangian with respect to the velocities is singular or, equiv- 
alently, when the Legendre transformation from velocity space to phase space is not 
invertible). Constraints naturally appear in the formalism as a consequence of the 
equations of motion -if we are in velocity space- or due to the fact that the Legendre 
transformation is not invertible -if we are in phase space. In phase space, the exis- 
tence of constraints is compulsory, in velocity space not. The confusion can arise when 
one considers regular (non-singular) Dynamical Systems which are deliberately con- 
strained to describe the motions in an ad hoc given surface. This surface can be either 
in configuration space (holonomic constraints) or in tangent space or phase space. In 
the holonomic case, one usually deals with these systems with the theory of Lagrange 
multipliers[p [Q] [Q] which is physically based on D'Alembert's principle. Let us 
notice that in general the presence of these ad hoc constraints ammounts to a change 
of the original equations of motion. 

Dirac ||5|[^, in his pioneering work, studied CDS in order to get a Hamiltonian for- 
mulation of gauge theories, including General Relativity (GR). As a general covariant 
theory, GR contains some gauge transformations (i.e.: symmetries that depend upon 
arbitrary functions of space-time. In the case of GR they are the spacetime diffeomor- 
phisms) and it turns out that the price for having this kind of transformations in the 
formalism is that the Lagrangian must necessarily be singular. In fact, not only GR, 
but the most important quantum field theories, and also string theory, have also room 
for gauge transformations. This fact makes CDS a central issue in the modern study 

^The canonical formalism is built out of it by using the Legendre transformation. 
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of Dynamical Systems. 

The existence of gauge transformations means that there are unphysical degrees 
of freedom in the formahsm. This is also reflected in the fact that there is some 
arbitrariness in the dynamics, since to a given set of initial conditions there correspond 
several -actually, infinite- solutions of the equations of motion, which are related among 
themselves through gauge transformations. To get rid of these transformations, i.e., 
to quotient out the spurious gauge degrees of freedom, we must somewhat reduce the 
dimensionality of tangent space or phase space (depending upon we are working in 
Lagrangian or canonical formalism). One way of doing that, which proves convenient 
in varied circumstances, is by ad hoc introducing a new set constraints in order to 
eliminate the unphysical degrees of freedom. These constraints are called Gauge Fixing 
(GF) constraints and its role is twofold Q Q|: fixing the dynamics and setting the 
physically inequivalent initial conditions. Notice that now we are introducing ad hoc 
constraints into the formalism as it can be done for regular theories. The difference is 
that now we do not intend to modify the dynamics, but rather to fix it by selecting 
one specific dynamics from the -gauge related- family of possible dynamics described 
by the equations of motion. 

In this paper we will consider the GF procedure to fix the dynamics of a CDS in 
the case of an holonomic GF (i.e.: constraints defined in configuration space Q). This is 
obviously the simplest case and it allows for both Hamiltonian and Lagrangian analysis. 
We find it instructive to compare the role of these ad hoc GF constraints introduced 
in a CDS with the role -dictated by D'Alembert principle- of the ad hoc constraints 
introduced in a regular Dinamical System. The geometric version of D'Alembert prin- 
ciple is: The holonomic constraints define a reduced configuration space Q-ji which has 
a natural injective map to Q, i : Q-ji — > Q. i is naturally lifted to the tangent map 
between the tangent bundles, i' : TQ-ji — > TQ. Then, the new dynamics is defined in 
the velocity space TQ-ji through a Lagrangian L/j which is the image of L (the original 
Lagrangian) under the pullback i'* of i' . In plain words this means that we get just 
by substituting the constraints into L. 

This is standard theory for constrained regular Dinamical Systems, of course. But 
some questions arise when we try to use the same mechanism in order to fix the dy- 
namics (through GF holonomic constraints) of a CDS. Now, in the singular case, we 
are lacking of any physical principle -like D'Alembert's- to justify this procedure, but 
since it is available, it is worth to explore it. The question is: Can we proceed in the 

^this is only a part of the whole GF procedure. The second part, as we have just said, consists in 
fixing the initial conditions. 
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same lines as it is done in the regular case, i.e.: to produce a reduced Lagrangian which 
dictates what the dynamics shall be?. In other words: it is correct to substitute the 
constraints into the lagrangian in order to get the correct gauge fixed dynamics? Is 
there any loss of information of the original theory if we proceed this way?. Our answer 
will be that the procedure is correct but, since there really is a loss of information, it 
must be supplemented with the addition of a specific subset of the primary Lagrangian 
constraints of the original theory. 

Throghout the paper we will work mainly in tangent space but we will turn to 
the canonical formalism when we find it convenient. Some standard mathematical 
conditions for the Lagrangian are assumed, namely: the rank of the Hessian matrix 
is constant and -in Hamiltonian picture- no second class constraint can become first 
class through the stabilization algorithm. 

In section 2 we will briefly consider the reduced Lagrangian formulation of the 
system. This formulation is equivalent to the extended one when use is made of the 
Lagrange multipliers. We will show in particular that the Lagrange multipliers for a 
CDS GF constraints are combinations of the original primary Lagrangian constraints 
of the theory. 

The main result of this paper is (section 3) the following: Once the holonomic GF 
constraints to fix the dynamics are introduced, a) The reduced Lagrangian is singular 
if and only if the original theory has Hamiltonian second class constraints, and b) there 
is no gauge freedom for the reduced Lagrangian (i.e.: it has no first class Hamiltonian 
constraints). 

In section 4 we illustrate our results with some examples. Comments and conclu- 
sions are given in section 5. 

2 formalism in the reduced velocity space 

Let us first consider a time-independent Lagrangian L(q,q) in TQ, where q = {q^,i = 
1, • • • , n} are local coordinates for a point in a n-dimensional configuration space Q. At 
this moment L can be either regular or singular. 

Let us introduce a set of -independent- holonomic constraints f^iq) = 0, = 
1, - ■ ■ ,k < n. These constraints define a reduced configuration space Q-ji where we can 
define coordinates Q = {Q"',a = 1, ■ ■ ■ , N := n — k}. The injective map Q-jz Q is 
defined by some functions = q^{Q^) such that the rank of \dq^ ldQ"-\ is k (maximum 
rank), and the lifting of this map to the tangent structures allows to define the reduced 
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Lagrangian L/j as 

Lr{Q,Q) :=L(g(Q),^Q). 

Then, working in the second tangent bundles, it is easy to see that ^ 

dq'' 

[^«]« = W^ag^' 

where [L]i stand for the functional derivatives of L: 

r -, dL d , dL , 

Since df^/dq^ form a basis for the independent null vectors of dq^ /dQ""^ i.e., 

^^-0 (2 1) 

identically, we can conclude that there exists some quantities such that the equations 
[LR\a = are equivalent to [I/]j+A^(9/^/9g* = ^ , fu{q) = 0, which in turn are equivalent 
to the equations obtained from the variational principle for the extended lagrangian 
Le := L + A^/^, where A'^ are taken as new dynamical variables. 
So we have the well known result 

[Lr] = ^ [Le] = 0, 

which is independent on whether L is singular or not. Let us first consider the regular 
case. The evolution operator in TQ derived from L^; is 

where 

x„ = ^|-+„.(,y-'r^|-. 

Stability oi fy := X./^ = X.^/^ under X determines the Lagrange multipliers as 
where 9^^'^ is the inverse of 



^To alleviate the notation we sometimes do not write the puUbacks explicitely. Here for instance 
the puUback of the functions in T^Q to functions mT^Q-iz is understood. 
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Observe that, except for the case when is alreay stable under Xq, we end up 
with a dynamics which is different to the original one. 

Now for the singular case. The equations of motion obtained from L are: 

- W.jii^ = 0. 

Since Wij is singular, it possesses r null vectors 7* , giving up to r (independent or 
not) constraints 

= . 

It proves very convenient to use a basis for these null vectors which is provided from 
the knowledge of the r primary Hamiltonian constraints of the theory, (f?p. Actually 
one can take [^: 

7j = ^b.p). (2.2) 

where Pi{q,q) = dL/dq^. It is easily shown that there exists at least one M'-? and 7^ 
such that 

and therefore |10| 

q' = M^-'a, + rYp, 

where fj'' can be taken as arbitrary functions of t. 

The stabilization algorithm starts by demanding that time evolution preserve the 
constraints (Xijp. Sometimes new constraints are found; sometimes some of the fjP are 
determined; eventually the dynamics is described by a vector field that exists on, and 
is tangent to, the constraint surface in velocity space: 

d d d 
+ + 



the a' are determined from the equations of motion and the stabilization algorithm; rj^ 
(fi = 1, ■ ■ ■ ,pi) are arbitrary functions of time; and 
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ir^M = Kjr^l = =■■ \D^u\ + 0. (2.3) 



where 7^ are a subset of the null vectors of Wjj, corresponding to the first class primary 
constraints 0^ Q found in the Hamiltonian formalism |^ (which we take in number 
r\ < r). 

A holonomic GF for the dynamics will consist in the introduction of r\ independent 
functions /^((7) in such a way that the arbitrary functions rf" (fi = 1, • • • , ri) of the 
dynamics become determined by the requirement of stability of fu ■= X/jy = ^ofu 
under the action of X. To do so we need: 

(observe that the relation |{/!/, 0^}| 7^ shows that the GF constraints = also fix 
the dynamics in canonical formalism, since 0^ are the first class primary constraints in 
the Hamiltonian formalism.). 

With these GF constraints f^{q), the extended Lagrangian Le gives the equations 
of motion: 

Now compute A'^. Contraction of 7* with the first set of these equations gives 



(we have supressed some coordinate indices in the last expression) whereby we can get 
as 

A'^ = -(D-^)^'''a-f^.. (2.4) 

The noticeable fact is that now the Legendre multipliers A'^ are constraints. This 
is good news because it tells us that the introducction of the GF constraints has not 
modified the dynamics. Let us be more specific on this point. 

There is a splitting of the set of hamiltonian primary constraints as first and second 
class (we will use indices ^ for first class, and fi' for second class) constraints^ It will 
prove convenient later, since the first class constraints (/>^ satisfy \{fu, ipli}] / 0; to take 
a basis (pjj^, for the second class constraints such that 

Wl,fu} = 0; (2.5) 



*Here we refer to first class primary constraints as the subset of primary constraints that still satisfy 
the first class condition after we have run the stabilization algorithm to get all the constraints of the 
theory: primary, secondary, etc. If we only look at the primary level, the number of first class primary 
constraints can be greater. 

^The previous footnote applies here. 
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then 

7^'|^ = 0. (2.6) 

This spUtting is translated to the Lagrangian formalism through ( |2.2| ), so we 

end up with primary Lagrangian constraints ~ and 07^/ ~ 0. It is worth noticing 
that only the first set is involved in the determination of A'^. 
Now observe that, due to (2.4): 
df 

[L] + = , 07^ = ^ [L] = , 07^ = , ^ [L] = 

the last equality holding because the constraints a^^ = are consequence of [L] = 0. 
Therefore 

[L]=0, f^{q)=0^ 
W+A^^=0,a7^ = 0,^(<?)=0 
^ [Le] = , 07^ = ^ [Lr] = , 07^ = , (2.7) 

where 07^ is obviously understood with the pullback to TQn, i'*[a'^y). 

So we see that [L] = , f^{q) = =^ [Lr\ = but the converse is not true. In the 
next section we will get a perfect understanding of this fact. 



3 When is Lji regular? 

With the same notation as in the previous section, with L being a singular Lagrangian 
and Lr the reduced Lagrangian after an holonomic gauge fixing, we are going to prove 
the following 

Theorem 1. Lr is regular if and only if L has only first class (Hamiltonian) con- 
straints. 



3.1 Proof of theorem 1 

First observe that we can write 



Q^Lr dq^ dq- 



where Wij stands, as before, for d^L/d(fdqK We must check whether the Hessian 
matrix for Lr is regular or not. That is, look for the existence of solutions V"" of 
{d'^LR/dQ''dQ^)V^ = 0. Then, using (O): 
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and since df^/dq^ form a basis for the null vectors of dq^ /dQ^, there must exists rj^ 
such that 



''dQ- dq^ ' 



contraction with 7^, which are null vectors for Wij, gives: 

^'^ dqi - ^^/^ ' 
but since \Dy^\ ^ we conclude that ij^ = 0, which means 



M/«|;V" = 0. (3.2) 



Now, since the set of null vectors of Wij is 7^ , 7^, , from ( |3.2| ): 



for some 5'^ , 5'^ . Contraction of this last expression with df^/dq^, and use of (|2.1[) and 
Q, gives: 

^"7:?^ = <5"^.u = ^ 5^ = , 



where (^^) has also been used. Therefore 



for 5'^' arbitrary. This means that we will get as many independent null vectors V 
for d^Ln/dQ"^ dQ^ as indices run for v' , which is the number of second class primary 
Hamiltonian constraints. This proves the theorem, for if there are no second class 
primary constraints, there are no second class constraints at all. In such a case, there 
are no null vectors for d'^Lji/dQ"' dQ^ and Lr is regular. 

To obtain a basis for the null vectors V"" is is convenient to have a deeper look at 
the canonical formalism. 

3.2 canonical formalism 

Here we will only consider some results that prove interesting to us. There is a natural 
map 
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Momenta pi are mapped to momenta Pa of the reduced formalism according to Pa = 
Pidq^ / dQ"^ . Functions '4'{q,p) in T*Q are projectable to functions tpiQjP) in T*Qti if 
and only if 

HQiQ),p)=i'{Q,p^) (3.3) 



This projectabihty condition can be writen in a more famihar way. From ( p.SI ): 

dif; dip dq^ 
Wi " 'dPadQ^' 

Contraction with df^/dq^, and use of ( |2.1| ) gives 

dip dff, 
dpi dq^ 

This is the version of the projectabihty condition for ip we were looking for: 

=0 (3.4) 

So, according to (|2.5| ), the second class primary constraints (/)^/ are projectable to some 
functions 0^/. We will prove that 0^, are the primary Hamiltonian constraints coming 
from the reduced Lagrangian Lr. First we can see that they are in the right number 
because it coincides with the number of null vectors V we have found for the Hessian 
matrix of Lr. Thererfore we only have to check that the pullback of (pj^, to the tangent 
space TQ-ji is identically zero (this is the definition of primary constraints). In fact: 



identically. We have thus got a basis for the null vectors V as = d(j)j^,/dPa := 7^/- 
Now we can get the primary Lagrangian constraints for Lr. It is easy to see that 

OLr O^Lr _. , . d^qi , dq' _ 



then the primary Lagrangian constraints aR7^' for Lr are 



dQ^dQ^^dpi dp 
where obvious pullbacks to TQn are understood. 
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We now have a complete understanding of ( p.?] ): 

[L] = 0, f^{q) = ^ [Lr] = 0, = 0, (3.5) 

for the dynamics for Lr only provides with the primary constraints api'y^i , whereas the 
rest of primary constraints for L, 07^, are absent in the reduced formulation, and must 
be separately introduced in order to maintain equivalence with [L] = , f^{q) = 0. 
Now we are ready to prove: 

Theorem 2. The dynamics derived from Lr has no gauge freedom. 

This is equivalent to say that all Hamiltonian constraints for Lr are second class. 



3.3 Proof of theorem 2 

First we will get the canonical Hamiltonian for Lr. The dynamics determined by 
the GF = is described in the canonical formalism by a first class Hamiltonian 
Hpc which can be taken to satisfy {f^^Hpc} = (Since the stabilization of = 



determines Hpc)- According to (3^), Hpc is projectable to a function Hr := Hpc in 

We can prove that Hr is a canonical Hamiltonian corresponding to Lr. First 
consider the fact, very easy to verify, that the Lagrangian energy Er := QBLr/BQ — Lr 
satisfies Er := i'*{E), where E is tha Lagrangian energy for L. Now, defining FL : 
TQ — > T*Q as the Legendre map derived from L, and FLr : TQn — T*Q'ji as 
the Legendre map from Lr, then the following property holds: for any function ^ in 
T*Q projectable to ip in T*Qn, FL*j^i! = i'*FL*^. Then Er = i'*{E) = i'*FL*HFc = 
FL*j^Hfc: which proves that Hr := Hpc is a good Hamiltonian for Lr. 

Next, it is easy to see that for two projectable functions, , (see (|3.4| )), its Poisson 
Bracket (PB) is projectable and satisfies: 

{V^} = {^,|},j , (3.6) 



where { , }r stands for the PB in T*Qfi (The proof of ( |3.6| ) is immediate). 

Now we can realize that the stabilization algorithm for the second class primary 
constraints (pf^i in T*Q is exactly the same as for the primary constraints in T*Qti. 
Indeed, the time derivative of 4>j^, is: 

•^V' = {(t>li'^HR}R = {(l)j^,,HFc} =■■ $1' , 

and so on. But since the set of constraints (p'j^,, ■ ■ ■ (until the stabilization algorithm 
eventualy ends when no new constraints appear) is second class, we conclude, using 
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( |3.6| ) that so it must be for </'^/,0^/, • • •• Hence Theorem 2 has been proved. Notice 
that Hji actuahy is a first class Hamiltonian. 

The final picture for the reduction L — through the GF constraints f^ = 
has been established: In the reduced theory there is only place for the second class 
Hamiltonian constraints of the original theory. First class constraints have simply 
disappeared. Not only the primary ones tpj^ but also the secondary ip'^, etc. that 
arise through the application of the stabilization algorithm to ipj^. Consequently, the 
pullbacks of these first class constraints to velocity space, also disappear from the 
reduced velocity space. All the information carried by the first class structure has 
been erased. This applies in particular to the gauge generators, which are special 
combinations -with arbitrary functions and its time derivatives as coeficients- of the 
first class constraints. 



4 Examples 

4.1 Axial gauges in electromagnetism 

In the introduction we have mentioned the example of electromagnetism in n ■ A = 
gauges, n^^ = {nP,n),A'^ = (^'',A). After elimination of the Lagrangian becomes 
regular. Gauss law is missing. But there is a compatibility between this missing 
constraint and the reduced dynamics. Now the time evolution vector field is tangent 
to the Gauss law, and if our initial (for, say, time t=0) configuration of the field A 
satisfies VE = 0, then the constraint will be satisfied for all times. In QED we can 
use the reduced Lagrangian to write down the propagator for the photon in the gauge 
n ■ A = 0, but to settle the asymptotic initial and final states, we must require the 



fulfillment of Gauss law |11]. 

Let L = —\F^uF^^ be our starting point (Our metric is (+,-,-,-)). The only La- 
grangian constraint is 07 = 9o(9iA*) + /S.Aq = 0, the Gauss law. The time evolution 
operator, which only exists on, and is tangent to, the surface defined by the Gauss law, 
is: 

X = / d^yA^' , - / d^yidAdaA>') + ^AA^ — + / d!^y\(y,t)^ — , 

where A is the arbitrary function of the dynamics reflecting the existence of gauge 
freedom. The GF n ■ A = Q will fix A as 

A = -^(n-V(9^.4^) + A(n-A)). 
rv' 
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If the variable Aq is eliminated by using the GF constraint, we end up with the reduced 
evolution operator -always tangent to the Gauss law constraint- 

Now let us proceed the other way around, that is, instead of substituting the GF 
constraint n ■ A = into the equations on motion, let us plug n ■ A = into the 
Lagrangian L to get the reduced Lagrangian Lr[A, A]: 

with Foi = Ai — :^di{n ■ A) . This Lagrangian is regular and its time evolution operator 
is: 



J.M. Pons, Gauge Fixing into the Lagrangian 



14 



Xr 



1 



- I d'y {AA, + ^{dod, + ^A)(n • A) + {d, + ^doMA^)) 



(4.2) 



■ " ' ~ ' ■ " ' ■ " SAi (y ) 

As it was expected, there is no trace of the Gauss law. We can check, though, that 
the evolutionary vector field Xr, in (^]^) differs from the evolutionary vector field X in 



(4.1) by a term which is proportional to the Gauss law. Indeed: 



n-' 



1 



X - Xk = / -,(M8,A') + ^A(„ ■ A))^ 



This result explicitly exemplifies the relation (3.5). 



4.2 Pure Abelian Chern-Simons 

Consider the Abelian Chern-Simons 2+1 Lagrangian L 
volved). The primary hamiltonian constraints are 

.^'^ := vr'" - e^^PAp = 0. 

(fp is a first class constraint, whereas 0^, (j)^ are second class. There is a secondary first 
class constraint ijj := (9j(vr* + e^'^^ Aj) = which plays the role of the Gauss law in this 
case. 

Use of the GF constraint A*^ = to get the reduced Lagrangian gives: 

Lr = eO^^A^A,-, 

which is still singular. Its only primary constraints are second class: (j)^ := tt^ — A2 = 
0, := TT^ + Ai = 0. No secondary constraints arise and Gauss law has disappeared. 
Observe that we are verifying our theorem 2: the reduced Lagrangian Lr is still 
singular (because L has second class constraints) but has no room for gauge freedom. 



= ^e'^'^f'FpiyAp (no metric in- 



5 conclusions 

In this paper we have proved some results concerning the correctness of plugging a 
set of holonomic GF constraints into a singular Lagrangian L. Our result is: if the 
GF constraints are taken in such a way that properly fix the dynamics of the singular 
theory defined by L, then the reduced Lagrangian Lr only keeps the information of the 
second class Hamiltonian constraints of the original theory. All first class constraints 
have disappeared from the reduced formalism. Therefore, to maintain equivalence with 
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the dynamics defined by L we must add the primary first class constraints to the 
dynamics defined by L/j. It is remarkable that all the information carried by the first 
class structure has been erased. No only at the primary level, but at any level of the 
stabilization algorithm. This applies in particular to the gauge generators, which are 
made up of the first class constraints. 

This is the picture in Hamiltonian formalism. In Lagrangian formalism the con- 
straints erased through this procedure of plugging the GF constraints into the La- 
grangian are just the pullbacks to velocity space of the first class Hamiltonian con- 
straints of the original theory. 

Lorentz covariance of Yang-Mills type theories requires the existence of secondary 
Hamiltonian constraints (whose pullbacks to velocity space will be part of the pri- 
mary Lagrangian constraints). Also general covariance for theories containing more 
that scalar fields (like General Relativity) require these secondary constraints also. 
Therefore, in all these cases, the price for plugging the holonomic GF constraints into 
the Lagrangian will be the loss of the secondary first class constraints. The examples 
provided in the previous section not only show how this happens but also show the 
consistency of the reduced theory with the missing constraints. 

In our paper we have dealt with holonomic GF constraints only. The advantage of 
considering this case is that we have to our avail a reduced configuration space out of 
which both the reduced velocity and phase space are built. 

Let us finish with two more comments. First, the theory we have developped can be 
easily extended to cases where only a partial gauge fixing of the dynamics is performed. 
In such cases, there is still some gauge freedom left for Ln and only a part of the original 
primary first class constraints -and its descendants through the stabilization algorithm- 
disappear from the reduced formalism. Second, In the usual case (Yang-Mills, Einstein- 
Hilbert gravity, etc.) where the stabilization algorithm in velocity space has only one 
step, the surface defined by = (see Eq. (|3.5| )) is a constant of motion for Lr, but 
not a standard one, since 07^ = c 7^ 0) will not be, in general, a constant of motion. 
Our first example can be used to illustrate this aspect. 
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